Universal Behavior of Two-Dimensional '^He at Low Temperatures 
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On the example of two-dimensional (2D) '^He we demonstrate that the main universal features of 
its experimental temperature T - density x phase diagram [see M. Neumann, J. Nyeki, J. Saunders, 
Science 317, 1356 (2007)] look like those in the heavy-fermion metals. Our comprehensive theo- 
retical analysis of experimental situation in 2D ^He allows us to propose a simple expression for 
effective mass M*{T,x), describing all diverse experimental facts in 2D ^He in unified manner and 
demonstrating that the universal behavior of M* (T, x) coincides with that observed in HF metals. 
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One of the main purposes of condensed matter physics 
is to unveil the nature of the non-Fermi liquid (NFL) 
behavior in various strongly correlated Fermi-systems. 
These substances, such as high temperature supercon- 
ductors, heavy fermion (HF) metals, 2D electron and ^He 
systems are the objects of intensive studies leading to 
understanding of many-body effects and quantum phase 
transitions responsible for the NFL behavior. Heavy 
fermion metals provide iinportant examples of strongly 
correlated Fermi-systems In these compounds, be- 

ing f-electron alloys, a lattice of f-electron spins couples 
to the itinerant electronic system by s — / Kondo ex- 
change interaction. The properties of such systems are 
now hotly debated as there is a common wisdom that 
they are related to zero temperature quantum fluctu- 
ations, suppressing quasiparticles and giving rise to a 
quantum critical point (QCP), where the systems tran- 
sit to different magnetic ground states generating their 
specific NFL behavior [2, 13|] . On the other hand, it was 
shown that the electronic system of HF metals demon- 
strates the universal low-temperature behavior irrespec- 
tively of their magnetic ground state Therefore it 
is of crucial importance to check whether this behavior 
can be observed in 2D Fermi systems. Fortunately, the 
recent measurements on 2D '^He become available [1, Q. 
Their results are extremely significant as they allow to 
check the presence of the universal behavior in the sys- 
tem formed by '^He atoms which are essentially different 
from electrons. Namely, the neutral atoms of 2D ■^He 
are fermions with spin S* = 1/2 and they interact with 
each other by van-der-Waals forces with strong hardcore 
repulsion (due to electrostatic repulsion of protons) and 
a weakly attractive tail. The different character of inter- 
particle interaction along with the fact, that a mass of 
He atom is 3 orders of magnitude larger then that of an 
electron, makes '^He to have drastically different micro- 
scopic properties then that of 3D HF metals. Because of 
this difference nobody can be sure that the macroscopic 
physical properties of both above fermionic systems will 
be more or less similar to each other. 



The bulk hquid "^He is historically the first object, to 
which a Landau Fermi-liquid (LFL) theory had been ap- 
plied This substance, being intrinsically isotropic 
Fermi-liquid with negligible spin-orbit interaction is an 
ideal object to test the LFL theory. Recently 2D ^He 
sample has been fabricated and its thermodynamic prop- 
erties have been thoroughly investigated 0, [g] ■ Our anal- 
ysis of the experimental measurements has shown that 
the behavior of 2D ^He is pretty similar to that of 3D 
HF compounds with various ground state magnetic prop- 
erties. Because of van-der-Waals character of interpar- 
ticle interaction, '^He has a very important feature: the 
change of total density of '^He film drives it towards QCP 
at which the quasiparticle effective mass AI* diverges 

13] ■ This peculiarity permits to plot the experimental 
temperature-density phase diagram, which in turn can 
be directly compared with theoretical predictions. 

In this letter we show that despite of very differ- 
ent microscopic nature of 2D ■^He and 3D HF metals, 
their main universal features are the same, being dic- 
tated by a Landau quasiparticles paradigm. Namely, 
we demonstrate that the main universal features of '^He 
experimental T - x phase diagram look like those in 
HF metals and can be well captured utilizing our no- 
tion of fermion condensation quantum phase transition 
(FCQPT) 0, 0, [U, m, m based on the quasiparticles 
paradigm and thus deriving NFL properties of above sys- 
tems from modified LFL theory. The modification is that 
in contrast to the Landau quasiparticle effective mass, 
the ^He effective mass M*(T,x) becomes temperature 
and density dependent. We demonstrate that the uni- 
versal behavior of M*{T, x) coincides with that observed 
in HF metals. 

Let us consider HF liquid at T = characterized by 
the effective mass M* . Upon applying the well-known 
equation, we can relate M* to the bare electron mass M 
S [li] M*/M = 1/(1 - NoF\pF,PF)/i). Here Nq is 



the density of states of a free electron gas, is Fermi 
momentum, and F^(j>f,Pf) is the p-wave component of 
Landau interaction amplitude F. Since LFL theory im- 
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plies the number density in the form x = /Stt^ , we can 
rewrite the ampUtude as F^{pf,Pf) = F^{x)- When 
at some critical point x — Xc, F^{x) achieves certain 
threshold value, the denominator tends to zero so that 
the effective mass diverges at T = and the system un- 
dergoes the fermion condensation quantum phase tran- 
sition (FCQPT) The leading term of 
this divergence reads 

M*{x) _ B _ X 



(1) 



where M is the bare mass, Eq. ^ is valid in both 3D 
and 2D cases, while the values of factors A and B depend 
on dimensionality and inter-particle interaction [l3|. At 
X > Xc the fermion condensation takes place. Here we 
confine ourselves to the case x < Xc- 

When the system approaches FCQPT, the dependence 
M*(T, x) is governed by Landau equation 0, 
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where n(p, T, x) is the distribution function of quasipar- 
ticles. The approximate solution of this equation is of 
the form Hi 
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M*{T) 



M*{x) 



/3/(0)ln{l + exp(-l//3)} 



(3) 



where Ai > and A2 < are constants of order unity, 
(i = TM*{T)/pj, and /(O) - F^{xc). It follows from 
Eq. ^ that the effective mass M*{T) as a function 
of T and x reveals three different regimes at growing 
temperature. At the lowest temperatures we have LFL 
regime with M*(T) ~ M*{x) + aT"^ with a < since 
Ai > 0. This observation coincides with facts 0, 
The effective mass as a function of T decays up to a 
minimum and afterward grows, reaching its maximum 
Ml,,(T, x) at some temperature T^i^^iyx) with subsequent 



diminishing as T^"^/^ IJ, |l5|. Moreover, the closer is 
number density x to its threshold value Xc, the higher is 
the rate of the growth. The peak value Mj^ grows also, 
but the maximum temperature T,„ax lowers. Near this 
temperature the last "traces" of LFL regime disappear, 
manifesting themselves in the divergence of above low- 
temperature series and substantial growth of M* [x) . The 
temperature region beginning near above minimum and 
continuing up to T„^gj^(x) signifies the crossover between 
LFL regime with almost constant effective mass and NFL 
behavior, given by T~^/^ dependence. Thus the T^ax 
point can be regarded as crossover between LFL and NFL 
regimes. The latter regime sets up at T < Tmax, when 
M*(x) — > 00, giving rise to T~'^/'^ effective mass decay 

il,13- 
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FIG. 1: The phase diagram of 2D ^He. The part for z < 1 cor- 
responds to HP behavior divided to the LFL and NFL parts 
by the line Tniax(z) oc (1 — z)'^^'^ , where Tmax is the effective 
mass maximum temperature. The exponent 3/2 = 1.5 com- 
ing from Eq. (O is in good agreement with the experimental 
value 1.7 ± 0.1 g]- The dependence M*{z) cx (1 - z)~^ is 
shown by the dashed line. The regime for z > 1 consists of 
LFL piece (the shadowed region, beginning in the intervening 
phase 2 < 1 [Hi, which is due to the substrate inhomogeneities, 
see text) and NFL regime at higher temperatures. 



It turns out that M*(T, x) in the entire T and x range 
can be well approximated by a simple universal interpo- 
lating function similar to the case of the application of 
magnetic field [1, [l^, [l^ ■ The interpolation occurs be- 
tween LFL (M* cx T2) and NFL (M* cx T^^/S) regimes 
thus describing the above crossover. Introducing the di- 

we obtain the desired 



mensionless variable y 
expression 
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Here M^{y) is the normalized effective mass, ci and C2 
are parameters, obtained from the condition of best fit 
to experiment. Equation ([1]) shows that M^^ cx 1/(1 — z) 
and it follows from ^ that cx T'"^/^. As a result, 
we obtain 



cx (1 - z) 



3/2 



(5) 



We note that obtained results are in agreement with nu- 
merical calculations [3,111. M*{T) can be measured in 
experiments on strongly correlated Fermi-systems. For 
example, M*(T) cx C{T)/T cx S{T)/T cx Mo(T) cx x{T) 
where C(T) is the specific heat, S'(T) — entropy, Mo(r) 
— magnetization and x(T) — AC magnetic susceptibil- 
ity. If the measurements are performed at fixed x then, 
as it follows from Eq. (|4|), the effective mass reaches the 
maximum at T = Tmax- Upon normalizing both M*{T) 
by its peak value at each x and the temperature by Tmax, 
we see from Eq. ([4]) that all the curves merge into single 
one demonstrating a scaling behavior. 

In Fig. [U we show the phase diagram of 2D ■^Hc in the 
variables T - z (see Eq. ([T])). For the sake of comparison 
the plot of the effective mass versus z is shown by dashed 
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FIG. 2: The dependence of the eflective mass M* (z) on di- 
mensionless density 1 — z — 1 — x/xc- Experimental data 
from Ref. [1] are shown by circles and squares and those from 
Ref. are shown by triangles. The effective mass is fitted as 
M*{z)/M oc A + B/{1- z) (see Eq. (HJ), while the reciprocal 
one as M/M*{z) oc Aiz, where A,B and Ai are constants. 
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FIG. 3: The normalized effective mass AI^ as a function of 
the normalized temperature T /Tma.:x. at densities shown in the 
left down corner. The behavior M%r is extracted from experi- 



mental data for S{T)/T in 2D ^He [gl and 3D HE compounds 



such as CeRu2Si2 and CePdi- 
universal function (|4]). 



Rhx [H, [T3|, fitted by the 



line. The part of the diagram where z < 1 corresponds 
to HF behavior and consists of LFL and NFL parts, di- 
vided by the line Tniax(2) oc (1 — 2;)^/^. We pay attention 
here, that our exponent 3/2 — 1.5 is exact as compared 
to that from Ref. Q 1.7 ± 0.1. The agreement between 
theoretical and experimental exponents suggests that our 
FCQPT scenario takes place both in 2D ^He and in HF 
metals. The regime for z > 1 consists of low-temperature 
LFL piece, (shadowed region, beginning in the interven- 
ing phase 2 < 1 01) and NFL regime at higher tempera- 
tures. The former LFL piece is related to the peculiarities 
of substrate on which 2D '^He film is placed. Namely, it is 
related to weak substrate heterogeneity (steps and edges 
on its surface) so that Landau quasiparticles, being lo- 
calized (pinned) on it, give rise to LFL behavior [5|, |6|. 
The competition between thermal and pinning energies 
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EIG. 4: The dependence M]^{T /T-m&yi) at densities shown in 
the left down corner. The behavior M'^ is extracted from ex- 
perimental data for C{T)/T in 2D ^'He [1 and for the magne- 
tization Mo in 2D "^He The solid curve shows the universal 
function, see the caption to Eig. O 



returns the system back to FC state and hence restores 
the NFL behavior. 

In Fig. [21 we report the experimental values of effec- 
tive mass M*{z) obtained by the measurements on '^He 
monolayer Q. These measurements, in coincidence with 
those from Ref. show the divergence of the effec- 
tive mass ai X = Xc- To show, that our FCQPT ap- 
proach is able to describe the above data, we present the 
fit of M*{z) by the fractional expression M*{z)/M cx 
A + B/{1 — z) and the reciprocal effective mass by the 
linear fit M/M*{z) oc Aiz. We note here, that the linear 
fit has been used to describe the experimental data for 
bilayer ■^He Q and we use this function here for the sake 
of illustration. It is seen from Fig. [2] that the data in 
Ref. [?] (3 He bilayer) can be equally well approximated 
by both linear and fractional functions, while the data in 
Ref. [§] cannot. For instance, both fitting functions give 
for the critical density in bilayer Xc ~ 9.8 nm~^, while 
for monolayer Q these values are different - Xc = 5.56 
for linear fit and Xc = 5.15 for fractional fit. It is seen 
from Fig. [21 that linear fit is unable to properly describe 
the experiment Q at small \ — z (i.e. near x = Xc), 
while the fractional fit describes the experiment pretty 
good. This means that the more detailed measurements 
are necessary in the vicinity x = Xc- 

Now we apply the universal dependence (H]) to fit the 
experiment not only in 2D -^He but in 3D HF metals as 
well. M^ijj) extracted from the entropy measurements 
on the '^He film at different densities x < Xc smaller 
then the critical point Xc = 9.9 ±0.1 nm~^ is reported 
in Fig. [21 In the same figure, the data extracted from 
heat capacity of ferromagnct CePdo.2Rho.8 and AC 
magnetic susceptibility of paramagnet CeRu2Si2 17| are 
plotted for different magnetic fields. It is seen that the 
universal behavior of the effective mass given by Eq. ([4]) 
(solid curve in Fig. [HJ is in accord with experimental 
facts. All substances are located at x < Xr, where the 
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system progressively disrupts its LFL behavior at ele- 
vated temperatures. In that case the control parameter, 
driving the system towards its critical point Xc is merely 
a number density x. It is seen that the behavior of the 
effective mass M^(y), extracted from S{T)/T in 2D ^He 
(the entropy S{T) is reported in Fig. S8 A of Ref. 0) 
looks very much like that in 3D HF compounds. As we 
shall see from Fig. [5l the amplitude and positions of the 
maxima of magnetization Mo{T) and S{T)/T in 2D ^He 
follow nicely the interpolation formula ([4]). We conclude 
that Eq. (g]) allows us to reduce a 4D function describing 
the effective mass to a function of a single variable. In- 
deed, the effective mass depends on magnetic field, tem- 
perature, number density and the composition so that all 
these parameters can be merged in the single variable by 
means of interpolating function like Eq. ([U, see also Ref. 

a. 

The attempt to fit the available experimental data for 
C{T)/T in 3He Q by the universal function Mli(y) is 
reported below in Fig. [H Here, the data extracted from 
heat capacity C{T)/T for ^He monolayer [§] and magne- 
tization Mq for bilayer are reported. It is seen that 
the effective mass extracted from these thermodynamic 
quantities can be well described by the universal inter- 
polation formula We note the qualitative similarity 
between the double layer [B] and monolayer ^8] of ■^He 
seen from Fig. [H 



the maximal entropy (5'/T)inax versus the number den- 
sity X. They are extracted from the measurements of 
S{T)/T on 3He bilayer 0. The fact that both the left 
and right panels have the same behavior shows once more 
that there are indeed the quasiparticles, determining the 
thermodynamic behavior of 2D ^He (and also 3D HF 
compounds near their QCP. 

To conclude, we have described the diverse experimen- 
tal facts related to temperature and number density de- 
pendencies of thermodynamic characteristics of 2D "^He 
by the single universal function of one argument. The 
above universal behavior is also inherent to HF metals 
with different magnetic ground states. Also, for the first 
time, the amplitude and positions of the maxima of the 
magnetization Mo{T) and the entropy S{T)/T in 2D ^He 
as the functions of density have been analyzed on Fig. 5. 
We obtain the marvelous coincidence with experiment 
in the framework of our theory. Moreover, these data 
could be obtained for ^He only and thus they were inac- 
cessible for analysis in HF metals. This fact also shows 
the universality of our approach. Thus we have shown 
that bringing the experimental data collected on differ- 
ent strongly correlated Fermi-systems to the above form 
immediately reveals their universal scaling behavior. 

This work was supported in part by the RFBR No. 
08-02-00038, DOE and NSF No. DMR-0705328. 
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FIG. 5: Left panel, the peak temperatures Tmax and the 
peak values Mmax extracted from measurements of the mag- 
netization Mo in '^He [^. Right panel shows Tmax and 
the peak values (5/T)inax extracted from measurements of 
S{T)/T in ^He ^. We approximate T^ax oc (1 - z)^/^ and 
(S/T) 

max ax 

(X A/{1 - z). 

On the left panel of Fig. [SI we show the density de- 
pendence of Tmax, extracted from measurements of the 
magnetization Mq{T) on ^He bilayer The peak tem- 
perature is fitted by Eq. ([5]). At the same figure, we 
have also reported the maximal magnetization Mmax- It 
is seen that Mmax is well described by the expression 
Mmax oc (S'/T)max « (1 " z)-\ scc Eq. (P . The right 
panel of Fig. [5] reports the peak temperature Tmax and 
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